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Large-scale simulations on three-dimensional (3D) frustrated anisotropic XY model have been 
performed to study the nonequilibrium phase transitions of vortex matter in weak random pinning 
potential in layered superconductors. The first-order phase transition from the moving Bragg glass to 
the moving smectic is clarified, based on thermodynamic quantities. A washboard noise is observed 
in the moving Bragg glass in 3D simulations for the first time. It is found that the activation of the 
vortex loops play the dominant role in the dynamical melting at high drive. 
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Dynamical properties of current driven vortex 
mattcimm interacting with random pinning poten- 
tials in type II superconductors have attracted con- 
siderable att^tian- both experimentallya&ElCrt] and 
theoreticallyBlliJ'tllEj A better understanding of various 
nonequilbrium phases and phase transitions is essential 
for explaining the nonlinear current- voltage (I-V) charac- 
teristics observed in experiments, conducting on samples 
in an external magnetic fieldJalU which is known as the 
first evidence of the first-order vortex lattice melting. In 
addition, this problem is closely related to an important 
class of phenomena in condensed matter physics, such 
as dynamics of sliding charge-density waves (CDW) in 
quasi-one-dimensional conductors, I 1 \ Wigner crystals in 
a two-dimensional (2D) electron gasjl3 as well as driven 
interface in random medially 

Vortex matter shows the three regimes of creep, depin- 
ning, and flow in its transport characteristics, depending 
on the drive. In the flow regime, the periodicity in the 
direction transverse to motion leads to a novel phase: 
moving Bragg glass (BrG), based on an elastic trans- 
verse equatioii_Qf motion proposed by Giamarchi and Le 
Doussal (GL)E9 As the driven force is reduced, the ef- 
fective pinning strength becomes larger. It has been ar- 
gued that moving vortex matter may decay first into a 
moving smectidl2lllJ and then into a moving liquid. It 
can be further driven into a creeping BrG below the de- 
pinning threshold. Recently, these rj»oving phases have 
been observed britb exjppr.imentallytra and in numerical 
simmationsMm'BWEl 

A precise analytical description is very challenging, es- 
pecially in the transition regime, because-it is quite diffi- 
cult to deal with the topological defectala in the plastic 
flow. Even within the elastic approach, an analytical 
study of driven vortex matter is hampered by ctaaaamic 
nonlinearities such as Kardar-Parisi-Zhang termliZI that 
governs the vortex dynamics on large scale. The under- 
lying mechanism of the dynamical melting still remains 
open question. In recent years, several three-dimensional 
(3D) numerical simulations have been performed .EjcJ^Ea 



However, one common shortage is that the moving BrG 
holds out to arbitrary high drives, which is obviously con- 
tradictory to any real experiments. 

In the present Letter, we report new results of nonequi- 
librium simulations for vortex matter in anisotropic 3D 
systems with weak disorder. The dynamical melting from 
the moving BrG to the moving smectic is found to be first 
order. The washboard noise is observed in the voltage 
power spectra in the moving BrG. Our results suggest 
that the thermally activated vortices play a dominant 
role in the dynamical melting at high currents. 

Vortex matter in type II layered superconductors can 
be describeci-hy-the 3D anisotropic XY model on a simple 
cubic latticeE#3>Ea 

H = - Jjj cos(<fc - - Ajj), (1) 

m 

where 4n specifies the phase of the superconducting order 
parameter on site i, A^ = (27r/$o) f° A • dl with A the 
magnetic vector potential of a field B = V x A along the 
z axis. The random pinning potential is introduced in 
the coupling strength in the xy plane Jjj = Jq(1 + peij), 
where e^ 's are independently Gaussian distributed with 
zero mean, and unit variance, p represents the pinning 
strength.ES The coupling between the xy planes is J z — 
Jo/r 2 , (r is the anisotropy constant). This model is 
relevant to high-T c superconductors and artificially lay- 
ered superconductors. For the data presented below, we 
typically choose p = 0.04 which models weak pinning 
strength, 1/r 2 = 1/40, and the average number of vor- 
tex lines per plaquette / = Z 2 i?/$o = 1/20, where I is 
grid spacing in the ab plane. Our system size is L = 40 
for all directions. 

In order to study the transport properties, we incorpo- 
rate the Resistivity-Shunted- Junction dynamics in simu- 
lations. Realizing that the sum of supercurrents into site 
i can be expressed in terms of the derivative of Eq. (1) 
with respect to <pi, the dynamical equations for the 0's 
are readily derived by requiring the sum of currents into 
each site to vanish 
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where J ex t,i is the external current which vanishes ex- 
cept for the boundary sites. The rjij is the thermal noise 
current with zero mean and a correlator {f]ij{t)riij{t')) = 
2akBTS(t — t'). In the following, the units are taken of 
2e = Jq = h = a = 1. 

In the present simulation, a uniform external current 
I x along x -direction is fed into ihe system. The fluc- 
tuating twist boundary conditioned is applied in the xy 
plane to maintain the current, and the periodic bound- 
ary condition is employed in the z axis. In the xy plane, 
the supercurrent between sites i and j is now given by 
4% = Jij sm(6i - 8j - A i} - v i} • A), with A = (A x , A y ) 
the fluctuating twist variable and Oi = <j>i + • A. The 
new phase angle Oi is periodic in both x- and y-directions. 
Dynamics of A can be then written as 

A « = l3 E [4i+%']- J a, (3) 

<ij>a 

where a = x, y. The voltage drop is V = —LA X . 

The above equations can be solved efficiently by a 
pseudo-spectral algorithmEJ due to the periodicity of 
phase in all directions. The time stepping is done us- 
ing a second-order Runge-Kutta scheme with At = 0.05. 
Our runs are typically (4 ~ 8) x 10 7 time steps and the 
latter half time steps are for the measurements. For a 
given current, the simulations are started from high tem- 
peratures with random initial phase configurations, and 
vortex systems are gradually cooled down. The present 
results are based on one realization of disorder. For these 
parameters^ the equilibrium BrG melts to a liquid at 
T = 0.254E3 and the critical current I C (T = 0) is esti- 
mated to be around 0.095. We calculate the internal en- 
ergy e per flux line per layer, the helicity modulus along 
the z axis T z , and the vortex structure factor in the xy 
plane, together with the density of dislocations pd and 
vortex-antivortex pairs p av in the xy plane. The helicity 
modulus describes the superconducting phase coherence. 

As shown in Fig. 1, evident jump of e is observed 
around the melting temperature T m for I = 0.5, which 
clearly indicates a first-order dynamical melting. The 
helicity modulus T c is very sharply set up precisely at 
T m , consistent with the first-order nature. The vanish of 
T z above T m shows that the dynamical melting marks 
the loss of superconducting phase coherence along the z 
direction. 

To characterize the spatial order of moving phases, we 
show the structure factors in the vicinity of T m at var- 
ious currents in Fig. 2. At high currents / = 0.5 and 
0.3(3> ^c(0)), moving BrG's with six well-defined peaks 
are found just below T m , which shows the quasi-long- 
range translational correlations both perpendicular and 
parallel to the-flow direction, consistent with the predic- 
tions by GL.I13 By mapping out the trajectories of the 



moving vortices, we found a set of periodic coupled elas- 
tic channels, corresponding to peaks in the K y = axis. 
The peak in K y ^ is anisotropic. The half width along 
the flow direction is considerably smaller than that trans- 
verse to the flow, indicating that the positional correla- 
tion along the flow direction is much stronger than that 
transverse to the flow, which provides a solid base for 
the moving BrG theory,E2l where a long range order is 
assumed along the flow direction. 

One can see from Figs. 2 (a) and (b) that the mov- 
ing BrG melts sharply into a moving smectic just above 
T m where Bragg peaks only remain in-J^ a = axis, con- 
firming the proposal by Balents et alO The current de- 
pendent anisotropy of the smectic peak will be discussed 
later. Tho-jpoving smectic was unnoticed in a previous 
simulatior£3 using isotropic 3D XY model for higher vor- 
tex density with smaller system size. In addition, the 
orientation of vortex lattice was found to be not aligned 
with the direction of the motion,E3 inconsistent with the 
characterization of the moving BrG. 

Through extensive simulations in vortex flow regime, 
we universally find that the appearance of a moving 
smectic is simultaneously accompanied by the loss of the 
superconducting phase coherence along the z direction 
(T z = 0). It demonstrates that the vortices flow incoher- 
ently in different xy planes in a moving smectic. 

As the temperature further increases, the moving 
smectic continuously evolves into a moving liquid with a 
ringlike pattern in the structure factor (not shown here). 
We have not found any anomaly in e, demonstrating a 
continuous phase transition. In low current regime, e.g. 
/ = 0.05 (< 4(0)), it is found in Fig. 2(c) that the 
(creeping) BrG directly melts to a liquid without through 
an intermediate moving smectic, close to equilibrium. 

To capture insight of the underlying mechanism of the 
dynamical melting, we display temperature dependence 
of the density of dislocations pa in Fig. 3. Sharp jump of 
Pd right at T m demonstrates that the dynamical melting 
is mediated by the proliferation of dislocations, analogous 
to the equilibrium meltingpll. Dislocations can be created 
either by entanglement of the field-induced flux lines or 
by excitations of small vortex loops, the latter could be 
more important in the presence of external currents. For 
/ = 0.3, we find that vortex loops are only activated 
far above T m , which excludes the effect of vortex loops 
in the melting. Anisotropy of the smectic peak at this 
current exhibited in Fig. 2(b) is therefore quite similar 
to that in Ref. (19) without the mechanism to generate 
vortex loops. However, at / = 0.5, the thermally ac- 
tivated vortices play an essential role in the dynamical 
melting. The density of thermally excited vortices p av 
also shows a steep jump at T m . The smearing out of the 
smectic peak transverse to the flow at / = 0.5 shown in 
Fig. 2(a) is just originated from the nucleation of vor- 
tex loops. As the current increases to 0.69, "thermally" 
activated vortex-antivortex pairs remain down to zero 
temperature. Above this current, the superconducting 
coherence along the z direction vanish at any tempera- 
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ture. In this sense, we refer this current to a depairng 
critical current. It is just the thermally activated vortex 
loops that induce dislocations, which in turn destroy the 
superconducting moving BrG in very high currents. The 
perturbation of these excited vortex Jnop-s has not been 
considered in the analytical theories. E2lEilO In previous 
simulations using Langyejin idy.na.m.ics of a fixed number 
of interacting particles,BE3&3'BEa the effect of ther- 
mally activated vortices are absolutely excluded, so the 
moving BrG remains unreasonably in arbitrary high cur- 
rents. 

Repeating such simulations at various currents, we 
compose a dynamical phase diagram of vortex matter 
as a function of I and T in Fig. 4. The continuous phase 
transition from the moving smectic to the moving liquid 
is determined by criterion that the smectic peak height 
from the ring background is less than 0.05. T m (I) is al- 
ways lower than T m (0) and shows a nonmonotonic char- 
acteristics. This naturally bring about the possibility of 
the reordering of driven vortex matter with the increase 
of drives predicted theoretically. 

To show the temporal correlation in BrG phase, we 
calculate the power spectra of voltage noise s(a>) = 
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t dt\ . In Fig. 5(a), well-developed peaks 



appear at the washboard frequency and harmonics for 
/ = 0.3 at T = 0. The voltage noise spectrum for / = 0.5 
at T = is also presented in the inset of Fig. 5(a) for 
comparison. The voltage cross the sample is V x (t) oc nv y , 
where n is the vortex deasity and v y is the vortex ve- 
locity in the y direction.Ell The washboard frequency is 
given by ujq = v y /a, where a is the vortex spacing. Due 
to the same vortex density, we really find that the volt- 
age ratio (« 2.01) is exactly consistent with that of two 
washboard frequencies. As the temperature increases, as 
shown in Fig. 5(b), the peak at u>o remains in all moving 
BrG phases and suddenly disappears just above T m . The 
moving BrG phase is then characterized by the peak at 
luq. Interestingly, we find that the peak at uj = 5wo owing 
to the network discretization is gone in temperatures far 
below T m , which excludes the effect of the artificial lat- 
tice pinning in the phase transition. In the moving smec- 
tic, the short range of the position correlation along the 
flow direction totally destroys the temporal order of the 
moving vortices. To the best of our knowledge, the wash- 
board noise, observed recently in both conventional^ and 
high-T c superconductorsja has never been reported in np=. 
merical simulations except in 2D at zero temperature.Eil 
We believe this to be the first observation of the wash- 
board velocity modulation of driven vortex matter in 3D 
numerical simulations. 
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(1965). 

Fig. 1. Temperature dependence of internal energy 
e (cycles) and helicity modulus along the z axis T z 
(squares) at I = 0.5. 

Fig. 2. The vortex structure factors just above and 
below T m for (a) T = 0.5, (b) T = 0.3, and (c) T = 0.05. 

Fig. 3. Temperature dependence of the density of dis- 
locations pd and thermally activated vortex-antivortex 



pairs p av a,t various currents. 

Fig. 4. Dynamical phase diagram in temperature- 
current plane. Solid lines with open cycles: I s * order (dy- 
namical) melting from moving BrG to smcctic. Dashed 
lines: continuous phase transition from the moving smec- 
tic to liquid. 

Fig. 5. Voltage noise power spectra s(u) for / = 0.3 
at (a) T = 0, (b) T = 0.118, and for / = 0.5 at T = in 
the inset of (a). 
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Fig. 1 Chen et al. 



(a) 1=0.5 , T=0.162(left) and T=0.164 (right) 




(b) 1=0.3 , T=0.118(left) and T=0.120 (right) 




(c) 1=0.05 , T=0.236(left) and T=0.238 (right) 

Fig.2.Chen et al. 
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Fig. 3. Chen et al. 
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